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1.  Introduction 


Consider  a  monotone  relationship  of  the  expected  response  to 
the  dosage  in  bioassay.  A  critical  dosage  is  considered  here  to 
be  the  treatment  level  associated  with  an  expected  response  thre¬ 
shold.  The  problem  of  determining  critical  dosages  is  known  as 
the  calibration  or  the  inverse  regression  problem.  Many  studies 
can  be  found  in  the  literature  concerning  the  estimation  of  critical 

V  V 

dosages  for  linear  and  non-linear  regression  models.  (Krutchkoff 
(1967)).  Methods  of  estimation  of  critical  dosages  in  the  single 
sample  case  provide  either  point  estimates  or  confidence  intervals 
(Fieller  (1964)).  Among  the  methods  which  involve  multistage 
design  approach  we  find  the  stochastic  approximation  methods  and 
sequential  search  procedures  (Eiehhorn  and  Zacks  (1973)).  The 
present  study  considers  the  problem  of  determining  critical  dosages 
in  a  situation  where  a  large  number  of  regression  lines  are  avai¬ 
lable  from  many  related  assays,  but  each  regression  line  is  based 
on  a  small  number  of  observations.  Thus,  rather  than  estimating 
the  critical  dosage  for  each  assay  indivudually ,  the  inter-block 
information  provided  by  the  various  regression  lines  is  intilized 
to  increase  the  precision  of  the  estimates.  This  is  done  within 
a  Bayesian  framework.  More  specifically,  the  model  assumes  that 

in  each  assay  the  dosage  -  response  relationship  is  represented  by 

o 

a  linear  regression,  with  the  same  variance  o  around  all  regression 
lines  and  normal  distribution  of  errors.  The  threshold,  n,  is"  the 
same  for  all  assays.  If  and  k  =  are  the  true  . 

intercept  and  slope  of  the  k-th  regression  line,  the  k-th  critical 
dosage  is  =  (  n  -  °^)/^,  .  >  0  and  **  00  if  £  0.  ^The 
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prior  distribution  of  (a,g)  is  chosen  to  yield  a  negligible  prior 
marginal  probability  for  {g  <  0}.  The  Bayesian  model  assumes  that 
(a,» »  k  =  are  priorly  independent  and  identically  distri¬ 

buted  vectors  with  a  properly  chosed  prior  bivariate  normal  distribu- 

•V- 

based  on  the  observations  performed  in  the  k-th  assay.  On  the 
basis  of  this  LSE  a  posterior  bivariate  normal  distribution  is  deter¬ 
mined  for  (0^,6^)*  This  distribution  yields  a  predictive  normal  distri¬ 
bution,  given  (a^b^)  for  a  response  Y(5)  at  a  dosage  £.  (See  Aitchison 
and  Dunsmore  1975).  The  minimum  predictive  risk  estimator  of  the  critical 

A 

dosage  is  defined  as  the  dosage,  5^*  which  minimizes  the  predictive 

risk,  i.e.,  E{ (Y(C)-n)2 | (a^.b^) }.  The  suggested  estimator  of 

depends  on  the  LSE  ^ajc*'3jc)  and  on  the  parameters  of  the  prior 

bivariate  normal  distribution  of  (a.  .g,).  When  the  number  of 

k  k 

assays,  N,  is  large  an  empirical  Bayes  method  can  be  employed  for 
estimating  the  prior  parameters.  In  the  present  study  we  develop 
the  formulae  for  the  empirical  Bayes  estimation  of  the  critical 
dosages.  The  formula  obtained  resembles  somewhat  Stein-type  esti¬ 
mators  of  a  multivariate  mean  vector  (Zacks  (1971)).  The  procedure 
developed  in  the  present  study  is  applied  for  the  determination  of 
critical  concentrations  of  benzo-soluble  organic- extracts  in  air 
samples  taken  in  1963  and  1964  from  53  and  54  different  sites  in 
the  U.S.A.  These  organic  extracts  were  tested  for  their  toxicity 
in  a  series  of  photodynamic  bioassays  (Epstien  et.  al.  (1965)). 

The  toxicity  of  the  benzo-soluble  extracts  from  the  pollutants  is 


tion.  Let  (a^.b^)  be  the  least-square  estimators  (LSE)  of  (a^ 
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a  function  of  their  chemical  composition  and  concentration  in  the 
air.  The  chemical  composition  varies  (at  random)  within  a  site  and 
between  the  sites.  The  model  developed  in  the  present  study  was 
found  suitable  for  the  determination  of  critical  air  concentrations 
(dosages)  for  each  site.  These  critical  dosages  can  be  compared  with 
the  actual  concentrations  of  the  organic  extracts  in  the  samples.  When- 
every  an  air  sample  contains  organic  extracts  with  concentration  higher 
than  the  critical  dosage  evidence  exists  of  undesirable  toxicity  of 
the  air  pollution.  Similar  applications  can  also  be  performed  in 
other  areas  of  the  empirical  sciences. 

The  present  study  consists  of  five  sections.  In  section  2  we 
specify  the  statistical  model  and  the  Bayesian  framework.  The 
method  of  determining  critical  dosages  by  minimizing  the  predictive 
risk  is  provided  in  section  3.  Section  4  is  devoted  to  the  empirical 
Bayes  approach  when  the  number  of  assays,  N,  is  large.  Finally,  in 
section  5  we  present  the  application  to  the  analysis  of  the  photo¬ 
dynamic  bioassays,  for  the  determination  of  critical  concentration 
of  benzo-soluble  organic  extracts  in  air  samples. 

2.  The  Statistical  Model  and  the  Baysian  Framework. 

Consider  N  sets  of  biological  assays,  having  dose  response 

relationship  Yfa^)  “  +  sk\i  +  tki-  k  ’ 1 . "• 

i  «=  1 . n.  ,  where  e,  ,  is  a  random  variable  normally  distributed 

*  k  ki 

with  expectation  zero  and  varaince  o^.  The  regressors  (i*=l, . . ,n^) 


are  the  log-dosage  applied  at  the  kth  bioassay,  and  are 

the  linear  regression  parameters  based  on  the  observations  (y  . , 

Kl 

k  ■  1,...,N,  i  ■  l,...,n^.  The  model  assumes  that  the  variance  o2 
is  the  same  around  all  the  N  regression  lines.  Determine  the  common 
least  square  estimators  (LSE)  a^  and  b^  of  the  linear  regression 
parameters  and  the  variance  around  the  regression  line  s2.  Let 

tC 

s>2  denote  the  pooled  estimator  of  this  common  variance,  i.e. 

N  N 

(2*1)  sp  =  kIX  ^nk-2)sk^k=l^nk-2^* 

The  large  number  of  assays  considered  in  the  present  problem 

and  the  typically  small  error  variance,  a2  ,  provide  estimates  s2 

with  samll  standard  error.  'Accordingly,  we  develop  the  following 

Bayesian  model  under  the  assumption  that  a2  is  known  and  substitue 

s2  for  a2. 

P 


According  to  the  theory  of  least-square  estimation  in  normal 


models,  (a^.b^)’  is  a  random  vector  having  a  conditional  bivariate 

normal  distribution,  with  an  expectation  vector  (a  ,3,)'  and 

k  k 

covariance  matrix  E,  ,  where 

k 


(2.2) 


where 


designates  the  mean  log-dosage  at  the  k'th  assay  and  SDX^  = 

_ 

/_k?=i5Xki  "  xk'  *  The  Bflyesian  m°del  assumes  that  each  assay 

can  be  considered  as  a  random  sample  from  a  larger  population  of 

assays.  Accordingly,  we  assume  that  (a.  ,6,  )  follows  a  bivariate 

K  K 

normal  prior  distribution  with  prior  expectation  (a0>f$0)'  and 
prior  covariance  matrix,  T;  i.e.. 


(2.3) 


N 


Given  the  estimates  (a^,b^)’  and  the  posterior  distribution 

of  the  regression  parameters  a^-so  a  bivariate  normal 

distribution  (Zacks  1971,  Box  and  Tiao  1973)  with  expectation  vector 


/«*! 

ki 

■( 

\  +  T(i  +  T)_1  jk 


<2'4)  Kl  '  W+T<{k  +  T)' 

and  covariance  matrix 

(2.5)  vk  -  T  -  T’($  +  T)-1T. 

3.  Bayesian  Determination  of  Critical  Dosazes , 


The  critical  dosage  ,  for  the  k-th  bioassay,  is  defined 

as  the  value  of  x  for  which  the  expected  response  is  rj,  i.e.. 


,  k  a  1, • • • ,N  • 


(3.1)  5,  =  n  ■  ■  ^  ,  k  =  1 . N  . 

We  comment  here  that  in  practical  applications  of  the  model  we  assume 


that  all  >0.  The  Bayesian  framework  assumes  a  normal 


!U 

°i  ^  d 


II  Si 


-HI  £  « 


f 
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marginal  prior  (posterior)  distribution  of  g^,  which  extends  over 
negative  values  too.  This  problem  is  overcome  in  applications  if  the 
error  variance  o2  is  relatively  small  and  the  dosages  in  the  bio¬ 
assays  are  well  designed.  In  such  cases  the  posterior  variance  of  g^ 
is  often  sufficiently  small  so  that  the  posterior  probability  of 
negative  g^  value  is  negligible. 


In  classical  statistical  analysis,  is  frequently  estimated 

by  the  least-squares  inverse  regression  statistic 


(3.2) 


n  "  \ 
b. 


k  =  1, . . . ,N. 


Fieller's  theorem  (Fieller  (1944))  is  often,  applied  to  obtain  classical 
confidence  intervals  for  5^.  The  application  of  Fieller’s  theorm 
in  the  Bayesian  framework  is  not  compatible  with  the  definition 
(3.1)  due  to  the  interchange  in  the  role  of  parameters  and  statistics. 

We  consider  therefore  two  types  of  Bayes  point  estimators  of 

* 

One  is  obtained  by  substituting  in  (3.1)  the  posterior  estimates 
* 

and  g^  of  and  g^»  respectively.  Accordingly,  we  obtain  the 

(pseddo)  Bayes  estimator 


(3.3) 


n  -  a. 


k  —  l*p*  .*,14  • 


Notice  that  5^  is  not  a  Bayes  estimator,  since  it  does  not  minimize 
a  prior  (posterior)  risk.  We  introduce  a  proper  Bayes  estimator  of 
by  considering  the  value  of  x  which  minimizes  the  predictive 
risk  E( (Y(x)-n)2) .  More  specifically,  we  minimize  the  predictive 


expectation 
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(3.4)  Q  (x;!?k)  =  E{(Y(x)-n)2|  ?k)  , 

k  =  1,...,N.  E{ • |  ?k}  designates  the  expectation  with  respect  to  the 
predictive  distribution  of  Y(x)  in  the  k-th  bioassay.  In  the  present  case 

the  predictive  distribution  is  the  normal  distribution  with  mean  a,  +8,  x  and 

k  k 

variance  cr2  +  VR(x)  =  a2+  VRU  +  2xVkl2  +  x2VR22,  where  VR  , 
i,j  =  1,2  are  the  elements  of  the  posterior  covariance  matrix  VR. 

We  apply  here  the  loss  function  (Y(£R)-n)2  rather  than  (£R  -  £R)2 
since  the  posterior  expectation  of  (n-aR)/8R,  given  !?R  ,  does  not 
exist.  Thus 


(3.5)  Q  (*,  V  -  a’  +  VkU  +  2xvkl2  +  +  Co*+f£x-„>*  • 


The  minimization  of 
(3.6) 


(3.4)  with  respect  to 

*  .  *o 

5,.  - 


?k  = 


kl2  k 


1  +  w<2 


x  yields  the  estimator 

,  k  l,*tijN. 


A  (1-a)  level  predictive  interval  for  Y(|)  is  specified  by  the 
prediction  limits 


(3.7) 


P(£;ct)  =  a  + 


&  e± 


'l-a/2 


>/Q(b 


where  z,  is  the  l-a/2  fractile  of  the  standard  normal  distri- 

l-a/2 

bution.  It  is  not  difficult  to  show  that  the  three  different  esti- 

~  *  A 

mators  of  £R,  namely  £R,  CR  and  £R  are  consistant  ones,  as  the 
number  of  observations  nR  around  the  regression  lines  increase  to 
infinity  and  Sl))^  increase  to  infinity  too.  However,  questions  of 
consistency  and  asymptotic  efficiency  are  irrelevant  to  our  problem  since 
generally  we  are  concerned  with  cases  of  small  number  of  obervations  in  each 
bioassay.  For  this  reason  we  adopted  the  Bayesian  approach  to  compenseate  for  the 
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lack  of  accuracy  due  to  this  deficiency.  As  we  show  in  the  next 
section,  an  empirical  Bayes  approach  can  utilize  the  information 
obtained  from  the  large  number  of  different  assay  to  determine  an 
adequate  common  prior  distribution  for  the  analysis  of  the  indivi¬ 
dual  assays. 


An  Empirical  Bayes  Approach  for  Large  N. 


Generally  it  is  a  difficult  problem  to  determine  that  proper 
prior  parameters  for  each  regression  line.  However,  if  the  analysis 
consists  of  a  large  number  of  regression  lines  from  different  assays, 
and  if  it  is  plausible  to  assume  that  the  regression  parameters 
(a^,B^),  k  =  1,...,N,  constitute  a  random  sample  from  the  sample 
bivariate  (prior)  normal  distribution,  one  can  estimate  consistently 
the  prior  parameters.  More  specifically,  under  the  assumption  that 
the  (true)  regression  parameters  (cx^*^)*  ^  =  1>***>N,  are  inde¬ 
pendent  random  vectors  having  the  same  bivariate  normal  distribution. 


with  mean  (a0,8Q)  and  covariance  matrix  T  then 


(4.1) 


is  an  unbiased,  strongly  consistent  estimator  of  (Oq.Bq)1  having 

a  bivariate  normal  distribution  with  covariance  matrix  ^-(T  +  ^Y,  i,  ) . 

N  N^k=lTk 
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(4.2) 


C 


1 

N-l 


1 

N 


JN)(a,b) 


where  a'  =  (a, ,...,a„),  b'  =  (b,,...,b„),  I„  is  the  identity 
IN  INN 

matrix  of  order  N  and  is  an  N*N  matrix  of  l's.  Notice 

N 

that  the  total  covariance  matrix  T  +  —  I^-j^k  cornPose<^  the 

"within  variance”  component  anc*  t*ie  '  ^etween  variance" 

component  T.  Thus  as  in  the  common  components  of  variance  model 
(see  Graybill  (1976))  and  unbiased  estimator  of  T  is 

(4.3)  *-C-sJiA 


We  remark  that  if  the  design  matrices  of  all  the  N  assays  are 
the  same,  i.e.,  $k  =  $  for  all  k  =  1,...,N,  the  above  formulae 
simplify.  We  further  remark  that  (4.3)  may  be  negative  definite,  if 
the  "within  variance"  component  is  large  and  N  is  not  sufficiently 
large.  If  this  is  the  case,  one  has  to  apply  a  different  approach,  or 
use  biased  but  consistent  estimators  of  T.  Finally,  given  the 

A  /S  ^ 

estimates  (ot0,B0)  and  the  unbiased  estimator  T  one  can  determine 
an  estimate  of  V^,  namely 


(4.4) 


Vfc  =  T  -  f  ($k  +  T  ) 


-1 


kll 

Vkl2 

kl2 

\22 

and  substitue  its  elements  in  (3.6)  to  obtain  an  empirical  Bayes 
estimate  of  £  .  This  estimator  is 


-10- 


(4.5) 


|*  -  $  /ft 2 

^k  kl2/t5k 


1  +  Kiz'K 


where  and  3k  are  obtained  from  (2.4)  by  substituting  (cxQ ,  3q)  * 

A  A 

for  (a0,60)f  and  5^  ~  (T)-a^)/&^.  Notice  that  the  emperical 
BAyes  estimator  (4.5)  is  a  shrinkage  estimator  whenever 

A  A  k  A  A 

^  ^k^k22  is  the  case,  in  particular  when  Vj^  > 


In  the  foilwing  section  we  provide  a  large  scale  application 
of  the  empirical  Bayes  approach  described  above. 

5 .  An  Application  of  the  Model. 

• 

The  model  discussed  in  the  previous  sections  was  applied  to 
the  analysis  of  a  large  scale  photodynamic  bioassays,  performed  by 
Epstein  et.  al.  (1965),  for  the  purpose  of  evaluating  the  toxicity 
of  organic  extracts  from  atmospheric  pollutatants .  Air  samples 
were  collected  in  1963  and  1964  from  53  and  54  different  sites  in 
the  U.S.,  respectively.  The  benzo-soluble  organic  particles  were 
chemically  extracted  from  the  air  samples  and  the  atmospheric  con- 

centrations  [g/m  of  air]  were  recorded.  Proper  solutions  of  the 

•••  * 

organic  extracts  (O.E.)  were  tested  at  three  dilution  levels 
-4  -5  -6 

d  =  10  ,10  ,10  [g/mi.],  These  preparations  were  applied  in 

wells  including  30  cells  of  Paramecia  Caudatum.  The  measured  res¬ 
ponse,  called  th  LT90,  was  the  time  (in  minutes)  required  to  immo¬ 
bilise  90%  of  the  cells  under  ultra-violet  irradiation.  The  measu¬ 
rement  of  response  was  truncated  at  t^  *  90  minutes.  Response 
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values  over  90  minutes  are  therefore  unavailable,  neither  the  pro¬ 
potion  of  living  cells  at  the  time  of  truncation.  Four  replicas 
were  perfomed  at  each  dose.  Simultaneously,  the  LT90  was  measured 
on  a  standard  synthetic  benzo-a-pyrene  (BaP) .  41  complete  assays 

of  the  1963  data  and  54  of  the  1964  data  were  available  for  analysis. 

For  the  statistical  analysis  define  x  =  -log^  d  -  5.  The 
model  assumes  that  £n(LT90)  is  normally  distributed  with  mean 
a  +  $x  and  variance  o7 .  This  model  links  the  analysis  described 
later  to  the  theory  developed  in  the  previous  sections.  Bialik 
(1978)  verified  that  the  £n  LT90  versus  log-dose  regression  lines 
of  the  standard  preparations,  correspond  to  each  year  of  test  data, 
were  not  significatly  different.  Accordingly,  the  analysis  presented 
here  does  not  have  to  adjust  the  regression  line  of  each  site  for 
varying  experimental  conditions.  In  Table  1  we  present  the 
basic  response  statistics,  the  regression  statistic  and 
the  expected  LT90  corresponding  to  the  actual  concen- 
tation  in  the  air  for  sites  of  the  1963  samples.  The  re¬ 
gression  parameters  (a,  6)  of  different  sites  are  not 
expected  to  be  the  same  due  to  the  different  chemical 
composition  of  the  O.E..  Since  the  toxicity  of  the  organic 
pollutants  is  a  combination  of  their  chemical  composition  and  atmos¬ 
pheric  concentration,  Bialik  (1978)  introduced  a  measure  of  toxicity, 
AIRLT90,  which  takes  into  account  both  factors.  An  equivalent 
air-dosage,  AD,  is  defined  as  the  atmospheric  concentration  of 
the  O.E.  in  a  given  site  solved  in  1  mP  of  preparation.  Let 


[ 
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(5.1)  XAIR  =  -logi()(AD)  -  5. 

Then  the  corresponding  predicted  LT90  is  given  by 

(5.2)  AIRLT90  =  exp  (a  +  b.XAIR  +  o'2/2). 


The  XAIR  and  AIRLT90  of  the  various  sites  are  given  in  Table  1. 
Notice  that  all  the  XAIR  values  in  Table  1  fall  in  the  experimental 
domain  (-1  <  x  ^  1).  Accordingly,  the  AIRLT90  values  are  not  based 
on  extrapolation.  We  also  remark  that  intensive  photodynamic  acti¬ 
vity  is  associated  with  low  LT90  values. 


In  Table  2  we  present,  the  Bayes  estimator  (a,  ,8,  )  and  the 

fC  tc 

-  *  A 

corresponding  and  f estimators  of  the  inverse  regression 

parameters  corresponding  to  the  threshold  n  =  2.9.  This  value  of 

is  the  smallest  £n(AIRLT90)  in  Table  1.  The  Bayes  estimators  (a  ,B  ) 

k  k 

were  determined  according  to  the  empirical  Bayes  apt>roach,  described 
in  Section  4,  based  on  the  1963  and  1964  data.  The  LSE’s  (a  ,b  ) 

K.  K 

of  the  1963  data  yield  the  empirical  Bayes  estimates 


lao\ 

I 

3.8628^ 

(5.3)  °' 

ss  I 

Uoj 

63  ' 

^0. 9241; 

The  corresponding  covariance  matrix  is 

.1645 

(5.4)  (C) 

S3  ! 

t>  6  | 

.0354 

.0354 

.0723 


The  corresponding  covariance  matrix  for  the  1963  data  is 


(5.5) 


63 


f.0018  ,0017\ 

.0017  .0035, 


Thus,  the  empirical  Bayes  estimate  of  the  prior  covariance  matrix  T  is 


(5.6) 


'.1627  .0337  j 

,.0337  .0683/ 


We  computed  similar  estimates  based  on  the  1964  data  and  obtained 


(5.7) 


^0/  64 


3.5551 
,1.0673  , 


(T) 


64 


.2437  .0232 

,.0232  .0718, 


For  the  purpose  of  comparison  we  present  in  Table  2  for  each 
site  of  the  1963  samples  the  empirical  Bayes  estimates  corres¬ 
ponding  to  1963  and  1964. 

*  is 

Observe  that  the  values  of  (a.  ,B^)'  are  very  close  to 
(a^,b^)'  ±  (S.D.(a^),  S.D.(b^))'.  As  a  result  an(^  £.  have 

similar  values.  Also  .0001  £  |£  -  £  .0623  for  every 

*  is 

k  =  1,...,41.  However,  a  visible  effect  on  £  and  £.  is  demon¬ 
strated  by  changing  the  prior  distribution,  corresponding  to  changes 
in  the  experimental  conditions.  This  can  be  seen  when  the  empirical 
prior  distribution  based  on  the  1964  data  is  applied  to  the  1963 

if  is  is 

data.  Finally,  the  predicted  response  at  XAIR,  Y  (XAIR)  **  a  +3  XAIR, 
can  be  compared  with  the  y-fractile  of  the  predictive  distribu- 

A 

tion  at  namely 
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N  is  the  number  of  observations  around  each  line. 

Expected  LT90  =  exp{A  +  BX  +  o2/2}  . 

A  and  3  are  the  least  squares  estimates,  a  is  the  standard  deviation  around  the  regression  lines. 
AIRLT90  =  exp {A  +  BXAIR  +  c2/2},  XAIR  is  given  Table  2 


Table  2:  The  Actual  and  the  Critical  Atmospheric  Concentrations  of  Organic  Extracts  for  1963  Data 


o-or>  —  —  o(Df~nr>-oN(Dor>o-*Nrt(M<*o«rn  —  r)-«>»-ona)<DOoo-ONor>ntf> 
ncwhino*  —  co-ooodidcocd  —  no'-orgn^-r'nrjtsnf^owo- 
-tooDnoNon-<ontvO''Or>(M««--on(M»^icao''Hn-oniinn^)<rooo  —  nr^r^riN<» 
ft'OonnN  -a. a)  ffloon-oO'foirMDno(f)'OC'-o»-<-Nntn»NONH^MM»N 


—  -  99~00000N0--<-'000--00000->«-0--iN--.-<--00-- 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  t 


I  I  I  I  I  I  I 


•0C'D'0»-n(MD>f(BNffl<tM,)^^0"Cl0'NO-0lfi0l-0O-0-.ni|!NJ)*ifin0'-<>)’<f- 
r^noo»p)r'O0''-«  —  riioininococMrv^jcDcamo-CDnpiD-'O'roDNNinncD-ooDi^o 
•<otM,iONOn«o^M>'ainfii>t-<j'frKn'OniC''<in<ifiiinfi-o«'CDO-nNMnNn 
if)'<i*)(,)nf»-oa)CO(Dorj<o^’0(D(^oopifw<,)"<30'-o<t»<-'^pj<»-rj<»,c\ioo(>-<-o<j-n(N 


HHOOHOOOOODO^HHOOO^^OOOOO-tHHnriHNtlHHrtHOOHf* 

I  •  I  I  I  I  I  I  I  I  i  I  I  I  «<  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  !  I  t  I  I  I  I 


«t*,«i-oso)f>.o-£i*t'»'Oiftr*ir>*r'>ooi  —  no>  —  —  '-isno<i-o«»-nr'N  —  •oonoom*' 
o(D-«(3P)r)N  —  caNtnnrv-oorjN-oroo'O-p'CoriicotMnoocncD^jmtin-oof'riiin  — 
tvoo-niONon-<o»'Nc-<i<ir)4-ijMiir('Caio«'0-onin'OMnii)0'<nO"«nM) 
0'OinnnN-omotior)-oo»iD(DMDntiin'00'0»-«M,)^n<tf<o«>»<i<frtN 


--oo-oooooko---ooo-hooo-o--«o--n-«---ooh- 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  T  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 


(DOiDONnfioiiM)-<o^inMfi(Mo<Mii"Oin-onin(B'flN<trtr«<r»'OOHriM\i^- 

N«rnrjm«rcDT-tori(oinin^j-ONrvi>osN-mi>--<tan««-o-oin«ocv>OTO'0'Niri-oo 

*-oo-noi^on-<o«fr>'&--oui«i-in-o«r(M<r'OtDO^in'OftiininN<fa3o-'ni^cO'ONO 

in>oif)r)nh.-oiii(n(non-oo»CDCDNa3nnin'00-o«r*<-«(vp)»n«*-(MO(M-<'0'»'rcN 


h«OOhOOOOCNO-<hhOOO--<000-0-"hO-«I«-h--hOO-- 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  i  i  i  i  i  i  t  i  i 


<nmnin-oinooiCN»oriifiisnoO'(MMnMn'OrtnO'^<tniiD<f<nnrano-riiinN 
a>iri<»’>-o«rro«-<oMnir.  oivo>«tNtti'OND-D)i>'(MCD>oir)o<ir'-oaDtD'0«i-0‘riitnrso3o 
rtoo'oor'onnjo^-f'o-oniniooinnin'oaio  —  lo^rjinLON'rffio  —  nmcU'Orsin 
sn-onnr,)N'Oa>cD(Don-oo>taia)ScDr)(Mn-ooo<j'-H«Nr)«»-r,)<S’rjoni^>o«t’f»is 


H  H  O  O  ^ 
I  I  I  I  I 


OOOOONOrtHrtOOO^H000H0HHHOHrtNrtrtHHHOOHH 

T  i  i  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i  i  t  i  i  i  i  i  i  i 


«t*r*ONCDi>'O'0*r«r'0inr'-irwr,)O0i  —  so>  —  —  -  tsno«s-o>»or,'tM  —  -oontDO.* 
om-'CDnrjN  —  mrjinc'ir'-ooriif'.'Ocao-O'NcnCiicDCMnO'anm'OCDtfin.ooNriiin.* 
NOC'(voNOn«ooMm-(in-o<iM\iN’ai»o«<jninoMofflOHr)S'CnNfi 
in'OionnN-oiDtDCDon-ootcocDrvcoCDfjn'Oo-ov^-tNntnsi-niOMx'O'j-rjN 


^HOOHOOOOOnOHHHOOO^rtOOO^OHHHOHrtNHHHHHOO*<«4 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  t  I  I  I  I  i  I  I  I  I  I  I  I  I  I  I  I 


—  —  N(DniriiNtno«o«t')-<(\in«»'in«r«i-  —  r-'ttNoano'  —  coo 
—  Nr>-o-o«iD--ooc\i»,>ir>cD«r'0'»  —  •o--ot,jaji>'<MncDir>«i-i''Oa)o-  —  o  —  niO'NO«»'«t 
nO'0'OioN'fir>oino'inofi/N(VCDif)nfi/<or300f'>ooa30“inniO(hr)fvn<ro<n'fr' 
»<oooo-N«o-ooo-OHHOoi''000'<oooO'<0(iiif-ONon"nio 


O  I  O  O  I 


I  I  O  I  O  O  O  O  I 


OOI  IOOIOOOI  I  I  —  OOOIOIOOO 


<«NN0-N0'-0a)NC0'rcDw(M^or'rjfii(DfntMn«»--0'0'iftriiMN<ro-00'«rs0««»'-0 

'0«trjD-(D(D-0'rr'OM<t«ro'in<i(MO'OD'iriD-N<»->rii<t>orliChnir>r)r>^in<Tr>.«tincD 

-0MmH0'N(0tlMnOMMhO'0O(DN-00-ftlC-nim-tD((l-01l--0IMVNfiMO0-Oa) 


OOCO  —  O  —  O  —  OQ  —  OOO  —  OOOOOO  —  O  —  O  —  OOOOOO  —  O  —  0->0  —  O 


ajojrsO«nnin-«Nivoaj-or>D>-o<-n'<Nrt-ononNfv<fcoa)fiiD--N^mo'D'r30 

oco<rr>«-of"'C>'-<eMt'-0'>oir>a>NniCDroinr>.ON<toor)eocMOO-«»'Conor>(n-orii*o 

o-o«f-<r-o-o-o-i^)o-tM-<0'fMfv^»'>r)0'ONNa)tvcD(MODina)rv«frii«rr^«roo-noi» 


n»nn»nnn»nnn»n^nnn^nnnnnnn«nnn»»^»n*nnn«* 


•tor>>t-o«a)'ON'Or«ncDO«cDCM  —  P)N'*Nmo5o>c,)ONO'rto»>-or)(D«f'N«ir>o-n 
'0<tfMe-rocD-c<»NON«,<rcovfii'<0'0'0>r)0"-or)0'M«roruO‘»#inoo-«,»^^r)<co 
•0S0'-0xD‘0((D«rD-noN0'0‘O-0ajNr^o0‘0i0'0i^j  —  CDCo-oo-'OO'MNni-orJD'OCD 


oooo  —  o  —  o  —  oo  —  ooo  —  oooooo  —  o  —  o  —  oooooo  —  o  —  o  —  o  —  o 


cDoe'Oor)«n-Kr.-<a]io<t«»<«ooO'Or>inNOP)'Ory«rN-o»<o-ONO‘CD«N-9 

oro<r)r>-ONo-'NNO'-oinKNrua(Dor>-orsfl-oona)rii(T)0'<»-CDnor>r>-on«»-r> 

0“Q>r-<»'0'0<o^oo>M-*0'nirN^irjr)o-on(Na)r'a»«o)iocoo(«rN<N<<o»nM<' 


n*nn»nnn«nnn«ntnon»nnnnnfin*nnrnf«t'»n»nnnt» 


najr»‘fiioniop>-or)<rrsOtt‘CDCDn(Doofii-oo-om-<o>io-ON^np>*»«Nnr>r> 

■ors(ho.«D>niCD«»-<>-rtONO'0-oocDr^N^)0'ND'n(^}-a)a»-oo>-oo>MNfii>ONO‘oB 


■ors(ho^o>niCD«fD-rto>>»(Mhooa)NN.oo'ND-n(^)-a)a»-oo>-oo>MNfii'Oh>0‘OB 
bbdb-d-b-bd-dbb-bbbdbb-b-b-bbbbbb-b  —  b-b-d 


tDr»r>.<i<'nn(,)^Nr>ocD-o<0‘0«roor»«fflcx--»ni<-r»io-o-n-oc5'00‘N'00 
om«rmr)<ir\0'<Mso>oinQjcD«cocon5-ON^'Oonrartno»a)noo«nn(«»-o 
o>  o  *  —  ♦•0-03i«-0D>N«0'Nrs^»>n0'O(MN(nscDMCDiocDn«fM<«'r^^ - ~ 


tonocnnn^in 

cn<rr'*-oonrt<r 


ntnn«noo*nnr)fn»nnn»nnnnnon<nnn'K«tn»n*n*t 


nonooaoivonoNoooaiooiiiooocDOOooiONooaiaitDOiD^onoo 


«n»fi'OM»»o-nin*cMj»'ONnoHNnfMifOM»i>o-Nn*«cMh 

HHrtM««MM>iNnnnnnnnn44<t«wnnA«(icnAA 


I 


Q) 

05 

cd 

• 

PQ 

CM 

rH 

II 

cd 

u 

d 

•H 

u 

O 

•H 

4J 

CX 

B 

Q) 

a> 

• 

> 

w 

•H 

<u 

• 

V 

JC 

o 

cd 

4J 

t— I 

4-4 

Q) 

d 

o 

M 

o 

d 

/->» 

*0 

o 

(0 

0) 

•H 

Q) 

CO 

4-> 

00 

cd 

cd 

cd 

rd 

U 

CO 

4-1 

o 

/~N 

d 

■u 

CQ 

a> 

•» 

u 

8 

d 

CO 

o 

d 

o 

o 

U-i 

•H 

O 

M 

4-» 

•H 

cd 

CO 

cd 

u 

<D 

u 

4J 

rH 

d 

cd 

cd 

0) 

B 

3 

a 

•H 

4-» 

d 

4J 

a 

o 

CO 

cd 

o 

0) 

0) 

o 

CO 

•ri 

Q) 

4J 

U 

CD 

cd 

4-4 

PQ 

O 

cx 

CO 

0) 

u 

O 

JC 

d 

6 

u 

a> 

4J 

rH 

cd 

Q) 

cd 

U 

> 

rH 

cd 

iH 

cd 

d 

o 

i  s 

cr 

a> 

•H 

4J 

PQ 

•H 

Q) 

U 

•o 

• 

00 

o 

d 

Sf 

cd 

cd 

\D 

CO 

01 

05 

o 

-d 

CO 

>  kD 

rH 

*d 

4J 

PQ 

no 

a 

<v 

d 

JC 

u 

cd 

4J 

cd 

:  3 

CO 

CO 

<  w* 

< 

vD 

05 

* 

•> 

*H 

Ctf  -5C 

CO 

M 

w* 

:  45 

4-4 

•» 

< 

O 

>< 

4  mj* 

/"N 

✓-s 

rH 

CM 

CO 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fHTi.n  P.t.  Entered) _ 

REPORT  DOCUMENTATION  PAGE  befoIeVompletogTorm 

L  REPORT  NUMBER  |2.  GOVT  ACCESSION  NO.  3-  RECIPIENT'S  CATALOG  NUMBER 

_ Technical  Report  No.  38^ _ _ _ 

4.  TITLE  (and  Subtitle)  S.  TYPE  OP  REPORT  4  PERIOD  COVERED 


REPORT  NUMBER 


EMPIRICAL  BAYES  ESTIMATION  OF  CRITICAL  DOSAGES 
HAVING  SMALLEST  PREDICTIVE  RISK. 


Technical  Report 

4.  PERFORMING  ORG.  REPORT  NUMBER 


17.  AUTHORS) 


Ora  Bialik  and  S.  Zacks 


e.  CONTRACT  OR  GRANT  NUMBERf*.) 

N  00014-75-C-0529  ^ 
PROJECT  NR  042-276 


t.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

4 

Department  of  Mathematics  &  Statistics 
Case  Western  Reserve  University 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

OFFICE  OF  NAVAL  RESEARCH 
ARLINGTON,  VIRGINIA  22217 

U.  MONITORING  AGENCY  NAME  A  AODRESSfU  dlllereru  from  Controlling  Office) 


12.  REPORT  DATE 

August  15,  1979 

13.  NUMBER  OF  PAGES 

15 

IS.  SECURITY  CLASS,  (ol  I  ft  It  report) 


114.  DISTRIBUTION  STATEMENT  (ol  thle  Report) 


UNCLASSIFIED 

IS*,  declassification/ downgrading 

SCHEDULE 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED, 


1  17.  DISTRIBUTION  STATEMENT  (ol  tho  ebetrect  entered  In  Block  20,  II  different  Irom  Repot t) 


lie.  supplementary  notes 


If.  KEY  WORDS  (Continue  on  reveree  elde  II  neceeeery  end  Identity  by  block  number ) 

Linear  regression,  Bayes  estimates,  empirical  Bayes,  critical  dosages, 
calibration,  predictive  risk,  photodynamic  bioassays.  . 

/  ,  1 

<<  I G  +  “ I  +'  .  . 

o>'f  ‘a 

VOkACs”  RACT  (Continue  on  reveree  elde  II  neceeeery  end  Identity  by  block  number) 

^An  empirical  Bayes  procedure  is  developed  for  the  estimation  of  critical 
dosages  in  the  linear  regression  case.  If  Y  =  a  -fvBx  e  is  the  basic  linear 
model,  the  critical  dosage  is  defined  as  £(n)  =  (n-a)/B,  for  3>0.  A  new  type 
of  Bayes  estimator  of  C(n)  i fi-  derived  under  the  criterion  of  minimizing  the  pre¬ 
dictive  risk  E<(«  +  The  empirical  Bayes  procedure  provides  con-  | 

sistant  estimators  of  the  ptfTor  parameters  when  a  large  number  of  independent 
repetitions  of  the  expcriiicnjt  is  available.  The  methodology  is  developed  toy 


FORM 
I  JAN  73 


lDITION  OF  I  NOV  1.4  15  OBSOLETE 

S/N  0102*  (.r-  014-  4601 


tCCuRITY  CLASSIFICATION  OF  THIS  PACK.  (W»,.n  O.l.  Knl 


SECURITY  CLASS.  HCATION  OF  THIS  PACE  OOimn  Dmtm  Enttrtd) 


analyze  a  large  set  of  photodynamic  bioassays,  for  the  determination  of  critical 
air  concentrations  of  benzo-soluble  organic  extracts. 


-Uh 


S/N  0.02-  LK.  (ill-  6601 


SECURITY  CLASSIFICATION  Of  IMS 


